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ABSTRACT.  A  new  method  is  presented  for  calculating 
the  hydrodynamic  coefficients  of  axisymmetric  bodies. 
The  method  utilizes  an  extended  version  of  the  Douglas 
potential  flow  computer  program  for  axisymmetric 
bodies  to  calculate  the  potential  and  velocity  distribu¬ 
tion  along  the  body  required  for  the  determination  of  the 
forces  and  moment  acting  on  the  body.  Based  on  a  com¬ 
parison  with  experiment  limited  to  one  body  shape*  the 
method  yields  values  of  the  hydrodynamic  coefficients 
for  blunt-  or  flat-based  bodies  which  are  adequate  for 
most  engineering  work. 
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This  report  presents  a  new  method,  using  the  Douglas  potential 
flow  computer  program,  for  calculating  the  hydrodynamic  coefficients 
of  axisymmetric  bodies.  The  development  given  here  represents  the 
first  phase  of  this  work,  which  is  restricted  to  blunt-  or  flat-based 
bodies.  The  next  phase,  involving  an  extension  to  streamlined  bodies, 
will  be  presented  in  a  future  report. 
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NOMENCLATURE 

Equation  of  body  surface 
Axial  force  on  body 
Transverse  force  on  body 

Functions  describing  variation  of  vs  along  body 
Functions  describing  variation  of  <f>  along  body 
Unit  vectors  in  the  xf  y,  z  directions,  respectively 
Body  length 
Moment  on  body 

Nondimensional  derivative  of  moment  with  respect  to  w 

Nondimensional  derivative  of  moment  with  respect  to  q 

Nondimensional  derivative  of  moment  with  respect  to  w 

Nondimensional  derivative  of  moment  with  respect  to  q 

Unit  vector  normal  to  body  surface 

Fluid  pressure 

Angular  velocity  of  body 

Body  radius 

Radial  coordinate 

Body  fixed  orthogonal  curvilinear  coordinates 
Time 

Axial  velocity  of  body 
Fluid  velocity 
Velocity  of  body  surface 

Fluid  velocity  components  in  the  s,  n,  0  directions,  respectively 
Transverse  velocity  of  body 
Inertial  rectangular  coordinates 
X  coordinate  of  body  CG 

Nondimensional  derivative  of  axial  force  with  respect  to  u 
Body  fixed  rectangular  coordinates 
Z  coordinate  of  body  CG 
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Zj,  Nondimensional  derivative  of  transverse  force  with  re¬ 
spect  to 

Nondimensional  derivative  of  transverse  force  with  re¬ 
spect  to  q 

Z'w  Nondimensional  derivative  of  transverse  force  with  re¬ 
spect  to  w 

Z^  Nondimensional  derivative  of  transverse  force  with  re¬ 
spect  to  q 

P  Angle  that  a  body- surface  tangent  makes  with  the  axis 
y  Angle  the  body  axis  makes  with  respect  to  the  horizontal 

p  Mass  density  of  fluid 

<J>  Potential  generated  by  body 

SUBSCRIPTS 

u  Due  to  axial  velocity 

w  Due  to  transverse  velocity 

q  Due  to  angular  velocity 

s  In  direction  of  increasing  s 

n  In  direction  of  increasing  n 

G  In  direction  of  increasing  6 


v 
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INTRODUCTION 

Traditionally,  there  have  been  only  two  methods  for  estimating  the 
hydrodynamic  coefficients  of  axisymmetric  bodies.  For  bodies  having 
relatively  large  length-to-diameter  ratios  and  blunt  or  truncated  bases, 
slender  body  theory  can  be  used.  For  bodies  having  a  streamlined 
afterbody,  the  semiempirical  methods  of  Ref.  1  and  2  can  be  used,  pro* 
vided  the  shape  does  not  deviate  too  much  from  that  of  a  typical  torpedo. 
However,  for  blunt-based  bodies  having  small  length-to-diameter 
ratios,  or  for  streamlined  bodies  having  shapes  that  fall  well  outside 
the  empirical  data  of  Ref.  1  and  2,  there  exists  no  satisfactory  method 
for  estimating  hydrodynamic  coefficients.  To  fill  this  void  and,  hope¬ 
fully,  to  improve  upon  the  two  methods  mentioned  above,  the  work  pre¬ 
sented  here  was  undertaken. 

For  bodies  moving  with  a  large  axial  velocity,  small  angle  of  attack, 
and  small  pitch  rate,  so  that  no  large  amount  of  vortex  shedding  occurs, 
it  is  reasonable  to  look  upon  the  flow  around  the  body  as  being  made  up 
of  two  regions:  (1)  a  region  near  the  body  where  viscosity  and  turbu¬ 
lence  play  an  important  part,  i.  e.  ,  the  boundary  layer;  and  (2)  a  region 
outside  the  boundary  layer  where  the  flow  is  strictly  potential.  Bodies 
having  blunt  or  truncated  bases  normally  do  not  develop  thick  boundary 
layers;  hence  they  may  be  treated  as  though  the  flow  around  them  were 
purely  potential.  Of  course,  the  fact  that  the  flow  separates  from  the 
body  at  the  base  must  be  properly  accounted  for  in  the  treatment. 

Bodies  having  streamlined  afterbodies  normally  develop  thick  boundary 
layers  there,  so  that  any  treatment  of  them  must  analyze  both  the 
boundary  layer  and  potential  flow.  Although  the  ultimate  purpose  of 
this  work  is  to  facilitate  the  handling  of  the  case  where  a  thick  bound¬ 
ary  layer  is  present,  the  analysis  presented  in  this  report  is  restricted 
to  blunt-based  bodies  where  the  boundary  layer  may  be  neglected.  It 
is  planned  to  combine  a  technique1  for  analyzing  the  boundary  layer 
with  an  analysis  similar  to  that  given  here  to  handle  the  case  of  a 
streamlined  afterbody.  This  work  will  be  presented  in  a  future  report. 

In  seeking  a  method  to  calculate  the  potential  flow  about  an  axi¬ 
symmetric  body  of  arbitrary  shape,  it  soon  became  apparent  that  tht? 
method  developed  at  the  Douglas  Aircraft  Co.  (Ref.  4  and  5)  had  no 
close  rival.  Not  only  was  the  method  essentially  exact,  but  a  com¬ 
puter  program  was  available.  It  was  therefore  decided  to  build  the 
method  for  calculating  hydrodynami.  .  coefficients  around  the  Douglas 
potential  flow  solution. 


The  fir*  phase  of  this  boundary  layer  work  is  presented  in  Ref.  3. 
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EXTENSION  OF  THE  DOUGLAS  PROGRAM 

When  the  idea  of  using  the  Douglas  potential  flow  computer  program 
to  calculate  the  hydrodynamic  coefficients  of  axisymmetric  bodies  was 
conceived,  it  was  possible  to  use  the  program  to  obtain  values  of  the 
fluid  velocity  at  the  body  surface  for  a  stream  flow  parallel  to  the  axis 
and  a  uniform  crossflow.  Since  this  was  not  enough  information  to 
calculate  a  complete  set  of  coefficients,  a  contract  was  let  to  the 
Douglas  Aircraft  Co.  to  extend  its  theoretical  work  and  include  in  the 
computer  program  the  following  additional  capabilities. 

1.  Solution  for  nonuniform  crossflow 

2.  Calculation  of  the  perturbation  potential  at  the  body  surface  for 
axisymmetric  flows  and  crossflows 

3.  Calculation  of  the  perturbation  velocity  at  the  body  surface  for 
axisymmetric  flows  and  crossflows 

With  these  additions  to  the  program,  all  quantities  required  for  the  de¬ 
termination  of  a  complete  set  of  coefficients  were  available.  A  dis¬ 
cussion  of  theoretical  development  work  and  extentions  to  the  computer 
program  is  given  in  Ref.  6. 

The  term  "nonuniform  crossflow"  used  above  is  admittedly  not  ex¬ 
plicit.  More  precisely,  this  condition  may  be  defined  as  any  flow 
where  the  potential  generated  by  the  body  satisfies  a  boundary  condition 
of  the  type 

- I  =  (v„)n=0  =  f(s>  cos  6 

3n  / n=0 

The  terms  "perturbation  potential"  and  "perturbation  velocity"  are 
the  potential  and  velocity  arising  from  the  presence  of  a  body  in  a  mov¬ 
ing  mass  of  fluid.  They  are,  of  course,  equal  to  the  potential  and  ve¬ 
locity  generated  by  a  body  moving  in  still  fluid,  provided  the  relative 
motion  of  the  body  and  the  fluid  at  infinity  is  equal. 


PHYSICAL  ARRANGEMENT 

The  physical  arrangement  considered  is  shown  in  Fig.  1.  An  axi¬ 
symmetric  body  is  in  planar  motion  in  unlimited  fluid  with  an  axial 
velocity  u,  a  transverse  velocity  w,  and  an  angular  velocity  q.  The 
plane  of  motion  is  the  X-Z  plane  of  the  inertial  X,  Y,  Z  coordinate 
system.  Two  coordinate  systems  are  fixed  in  the  body:  a  rectangular 
x,  y,  z  system  originating  at  the  center  of  gravity  of  the  body,  and  an 
orthogonal  curvilinear  s,  n,  0  system  with  the  origin  of  the  s  coordi¬ 
nate  at  the  body  nose.  An  angle  y  is  used  to  define  the  orientation  of 
the  body  axis  relative  to  the  X  axis,  and  coordinates  X0  and  ZQ  are 
used  to  locate  the  position  of  the  body  CG  in  the  X-Z  plane. 
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Body  cross  section 
viewed  from  tail 


FIG.  1.  Physical  Arrangement. 


BOUNDARY  CONDITIONS 

The  required  condition  at  the  body  surface  is  one  in  which  no  fluid 
penetrates  the  surface  and  no  gaps  occur  between  the  fluid  and  the  sur 
face.  Hence  the  fluid  velocity  component  normal  to  the  surface  must 
equal  the  velocity  of  the  surface  normal  to  itself.  The  following 
equation  states  this  requirement. 


( v  n )  n=0  -  V  B  •  N  ( 1 

Vb  is  the  velocity  of  the  body  surface  and  N  is  the  unit  normal  vector 
to  the  body  surface.  N  is  given  by  the  relationship 

9B  9B  9B 
-  i  + -  j  + - 


grad  B 
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where  B  (the  equation  of  the  body  surface)  =  f(x,  y,  z)  =  0.  Since  the 
body  is  axisymmetric,  B  may  be  expressed2  as 


b  =  Vy2  +  z2  -  R<*)  =  0 

from  which  there  is  obtained 


8B 

8B  8R 

dR 

3x 

8R  8x 

dx 

8B 

y 

8y 

^y2  +  z2 

8B 

z 

3z  \y2  + 


(3) 


Substituting  Eq.  3  into  Eq.  2  and  making  use  of  the  fact  (see  Fig.  1)  that 
at  the  body  surface  dR/dx  =  -tan(3,  y  =  R  sine,  and  z  =  R  cos  0,  there 
results 


N  =  sin(3i  +  sin  0  cos(3j  +  cos  6  cos  (3k  (4) 

Using  ^'ig.  1,  the  velocity  of  the  body  surface  may  be  written  directly. 
For  axial  motion  u:  Vg  =  ui  +  Oj  +  Ok 

For  transverse  motion  w:  VB  =  Oi  +  Oj  +  wk 

For  rotation  q:  VB  =  qR  cos  0i  +  Oj  -  qxk  (5) 

Substituting  Eq.  4  and  5  into  Eq.  1,  the  boundary  condition  equation  is 
obtained. 

For  axial  motion  u:  (vun  )n-0  -  u  sin  (3 

For  transverse  motion  w:  (vwn)n=o  =  w  cos  (3  cos  0 

For  rotation  q:  (vqn)n=0=  q(R  sin (3  -  x  cos  (3)  cos  0  (6) 

The  subscripts  u,  w,  and  q  are  introduced  here  to  mean  resulting 
from  axial  motion,  transverse  motion,  and  rotation,  respectively. 


^  Writing  the  expression  for  B  in  this  manner  implies  that  R  is  a  single-valued  function  of  x.  How¬ 
ever,  it  is  easily  shown  that  the  resulting  expression,  Eq.  4,  is  completely  general  and  applies  equally 
well  to  bodies  where  R  is  a  multivalued  function  of  x,  such  as  flat-nosed  bodies. 
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FLUID  POTENTIAL 

The  potential  generated  by  the  body  will  be  made  up  of  three  parts 
arising  from  the  three  separate  motions  of  the  body:  u,  w,  and  q 

<b  u  4>w  ) 

Since  the  potential  generated  by  an  cxisymmetric  body  moving  along 
its  axis  of  symmetry  is  naturally  axisymmetric,  the  most  general  ex¬ 
pression  for  <J>U  is 


4>u  =  uFu(s,  n) 


(8) 


In  Ref.  5,  Hess  shows  that  the  potential  generated  by  an  axisym¬ 
metric  body  satisfying  a  boundary  condition  of  the  type 

( vu)n=o  =  f(s)  cos  8 

has  the  simple  form 


<J>  =  f(  s,  n)  cos  6 

From  Eq.  6,  therefore,  it  is  apparent  that  the  most  general  expressions 
for  <J>W  and  4>q  are 


d> w  =  wFw(s,  n)  cos  0 

4>q  =  qFq(s,  n)  cos  0  (9) 


PRESSURE  AT  THE  BODY  SURFACE 
BERNOULLI  EQUATION 

Since  the  body  is  not  in  steady  motion  relative  to  any  possible 
choice  of  inertial  coordinate  system,  the  nonsteady  Bernoulli  equation 
must  be  used  to^determine  the  pressure  at  the  body  surface.  Assuming 
the  convention  V  =  grad<}>  where  V  is  the  fluid  velocity,  the  pressure  at 
the  body  surface  is  given  by 


P  n=0  “ 


,p) 

,  9t  /  n=0 


ip(v2)n=0 


(10) 
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EVALUATION  OF  (3<t>/at)n=0 

The  84>/3t  in  the  nonsteady  Bernoulli  equation  is  the  time  rate  of 
change  of  the  potential  at  a  point  fixed  in  inertial  space;  i.  e.  , 

a<t>rx,  Y,  Z,  t) 

at 

However,  <j>  is  known  as  a  function  of  s,  n,  0,  u,  w,  and  q  where 

s  =  f(x,  r) 
n  =  f(x,  r) 
x  •■=  f<x.  Y,  Z,  t) 
r  =  f<X,  Y,  Z,  t) 

0  =  f(X,  Y,  Z,  t) 
u  =  f(t) 
w  =  f(t) 

q  =  £(t) 

Therefore,  the  chain  rule  must  be  used.  Making  use  of  Eq.  7,  there  is 
obtained 


where  the  dot  over  a  symbol  indicates  the  derivative  with  respect  to 
time;  i.  e.  ,  u  =  du/dt  etc. 
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Equations  6,  8,  and  9  allow  the  following  relationships  to  be  written. 

3<t>„\  /8F u\ 

-  =<vuS)„=0=u  -  =  ufu<®> 

9  s  / n^O  \  ds  /mO 


8<t>J  laF,, 

=  (vwsk=0  =  W1 - I  COS  0  -  Wfw(3)  COS  S 


3s 


ds 


94»q\  /3Fq 

-  =  <vqs>n=0=  - I  COS  0  =  qf  (s)  cos  9 

9s  /  n=0  \  9s  / n=0 

;9<M 

[ -  =  (vun)n=0  =  U  sinP 

9n  n=0 


8<t>v 

9n  / n=0 


=  (Vwn)n-0=  W  COS  P  COS  0 


ld\ 

-  I  =  (v  )n=0  =  q(R  sinp  -  x  cos  (3)  cos  0 

\  9n  / n=0 


9<t\ 


=  0 


30  /  n=0 


f9cf) 


W 


30  / n^O 


=  R( v  )  =  - w( F  )  r,  sin  0  =  - wh  (s)  tun0 

'  w  0  'n=0  '  w  'n=U  w '  ' 


/d<bq\ 

I - I  =  R  ( v  )  »  =  -  q( 

\  an  /  q  0  'n=0 

\  90  /  n-0 


Fq)n=0  Sin0  =  _Clhq(S)  sin0 


9M 

l—  =  <Fu>n^  =  hu<S> 

9u  /  n^O 


dd> 


W 


3w  f  n-0 


=  (Fw)n-Ocos  8  =  hw(s)  cos  © 


8«>q 

-  =  (F  )n=0  cos  0  =  h  (s)  cos  0 

9q  /n=0  q  9 


(12) 
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From  Fig.  1,  it  is  seen  that  s  and  n  are  functions  of  x  and  r  through 
the  relations 


x  =  xnose  *  Jo'  cosPds  +  n  sinP 
r  =  J*  sin  (3d s  +  n  cos  (3 


(13) 

(14) 


where  (3  ~  f(s).  Taking  the  partial  derivative  of  Eq.  13  and  14  with  re¬ 
spect  to  x  and  solving  simultaneously  for  9n/8x  and  9s/9x,  and  taking 
the  partial  derivative  with  respect  to  r  and  solving  simultaneously  for 
8n/8r  and  9s/8r,  there  is  obtained 

/8n\ 

-  =  sin  (3 

\  8x  / n^O 


-  =  -cos  (3 

9x  /n-O 


(-) 

\  9  r  I n=0 


cos  (3 


9  r  ln=0 


sin  (3 


(15) 


From  Fig.  1.  the  time  rates  of  change  of  ( 1 )  the  coordinates  of  the 
body  center  of  gravity  and  (2)  the  angle  between  the  body  axis  and  the 
horizontal  are 


dx°  +  • 

- =  u  cosy  +  w  sin  y 

dt 

dZo 

-  =  -u  s  in  y  +  w  cosy 

dt 

dy 

—  =  q  (16) 

dt 

Also  from  Fig.  1.  the  coordinates  of  a  point  in  the  inertial  X,  Y,  Z  sys¬ 
tem  in  terms  of  the  coordinate  system  in  the  body  are 
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X  =  XQ  +  x  cosy  +  r  cos  0  sin  y 
Y  =  r  sinO 

Z  =  Z  0  -  x  siny  +  r  cos  0  cosy  (17) 

Taking  the  partial  derivative  of  Eq.  17  with  respect  to  time,  there  re¬ 
sults 


dX0  dy  9x  dy 

0  - - x  sin  y  -  +  cosy  -  +  r  cos  0  cosy  - 

dt  dt  9t  dt 

I  90  9r\ 

+  siny|-r  sin  0  -  +  cos  0  - I 

\  9t  dt  I 


90  9r 

0  =  r  cos  0  —  +  sin  0  - 

9t  9t 


dZo 

dy 

9x 

dy 

0  - - 

x  cos  y  -  - 

sin  y - 

r  cos  0  sin \  - 

dt 

dt 

9t 

dt 

I  06  9r\ 

+  cosy  -r  sinB  -  +  cos  0  - I  (18) 

\  9t  9t / 

Substituting  Eq.  16  into  Eq.  18  and  solving  simultaneously  for  9x/9t, 
9r/9t,  and  90/9t,  there  is  obtained 

/3x\ 

I - I  =  -u  -  qR  cos  0 

\  9t  / n=R 


ldr\ 

I - I  =  - (w  -  qx)  cos  0 

\  9t  frr R 


90  \ 

dt  I r=R 


sin  0 

(w  -  qx)  - 

R 


(19) 


Substituting  Eq.  12,  15,  and  19  into  Eq.  11,  the  expression  for  the 
( 0<t»  /  9 1)  n;=0  is  as  follows: 


9 


I 
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^  9<t> 

1 3t  , 


n^=0 


(ufu  +  wfw  cos  0  +  qfq  cos  9)(u  cos  (3  +  qR  cos  (3  cos  0 

-  w  sinp  cos  0  +  qx  sin  p  cos9)  +  (u  s  in  p  +  w  cos  (5  cos  0 

+  qR  sinp  cos  9-  qx  cos  p  cos0)(-u  sinp  -  qR  sin  p  cos  0 

-  w  cos  p  cos  0  +  qx  cos  p  cos  0)  +  (-wh^  sin  0  -  qhq  sin  0) 

(sin  0  sin  0  \ 

w - qx - I  +  uhu  +  wh  w  cos  0  +  qh  cos  0  (20) 

R  R  / 


EVALUATION  OF  (V2)nm0 

The  square  of  the  fluid  velocity  will  be  equal  to  the  sum  of  the 
squares  of  the  three  orthogonal  components  vs  ,  vn,  and  . 

<V2)n=0=  rKs>„=0  +  <vwS>»=0  +  <V»=o)2 

+  t(vu»>n=0  +  <  Vw„>„-0  +  (V’-J2 

+  t(vu,)n=0+  (vw,)n=0+  (vq  )„J2  (21) 

Substituting  Eq.  12  into  Eq.  21,  the  expression  for  (V2)n^Q  is 
(V2)  =  (uf  +  wf  cos  0  +  qf  cos  0)2 

n=0  u  w  ^  q  ' 

+  [u  sinp  +  w  cos  p  cos  0  +  q(R  sinp  -  x  cos  P)  cos  0]2 

sin  0  sin  0  \ 2 

W - qhq -  (22) 

.  R  R  / 


TRANSVERSE  FORCE,  AXIAL  FORCE, 
AND  MOMENT  ON  BODY 


Utilizing  Fig.  1,  the  expressions  for  the  transverse  force  FT ,  axial 
force  Fa*  and  moment  M  on  an  elemental  surface  distance  of  the  body, 
ds,  can  be  written  directly.  The  transverse  force  is  determined  by 
integrating,  around  the  circumference  of  the  body,  the  component  of 
the  pressure  force  in  the  FT  direction. 


dFT 

ds 


-R 


COS  P 


pn^0  cos  0  d9 


(23) 
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The  axial  force  is  determined  by  integrating,  around  the  body,  the 
component  of  the  pressure  force  in  the  FA  direction. 

dFA  .2r 

-  =  -R  sin  p  /Q  Pn<)dO  (24) 

ds 


The  moment  is  determined  by  integrating  around  the  body  the  contribu¬ 
tions  arising  from  (1)  the  product  of  the  pressure  force  transverse 
component  and  the  lever  arm  x,  and  (2)  the  product  of  the  axial  com¬ 
ponent  of  the  pressure  force  and  the  lever  arm  R  cos  9. 


dM  2  2w 

-  -  Rx  cos  p  fQw  pn4)cos  Gd9-R2  sin  fQW  pn=0cos  ©  dO 


contribution  from 
transverse  pressure 
force 


contribution  from 
axial  pressure 
force 


(25) 


Substituting  Eq.  20  and  22  into  Eq.  10  yields  the  expression  for  the  pres¬ 
sure  on  the  body  surface,  pn=0  •  Putting  this  expression  into  Eq.  23,  24, 
and  25,  carrying  out  the  indicated  multiplications,  and  integrating  with 
respect  to  0  around  the  body,  and  with  respect  to  s  along  the  body  from 
the  nose  (s  =  0)  to  the  tail  (s  =  sTAIL  ),  there  results 

Ft  =  puwrr  /JTAIL  R  cos  (3(- sin  (3  cos  p  -  fu  sinp  +  fw  cos  p  +  fufw  )ds 

+  puqTT  J*gTAIL  R  cos  p(-R  sin2  p  +  x  sinp  cos  p  +  Rfu  cos  p 

+  xfu  sinp  +  f q  cos  p  +  fufq)ds  +  pwn  f*TA1L  R  cos  phwds 

+  p qTr  jJTAIL  R  cos  ph^ds 

Fa  (apparent  mass)  =  pu2iT  /JTAIL  R  sinphuds 

M  =  puww  /qTA11  R(R  sinp  -  x  cos  P)(- sinp  cos  p  -  fu  sinp 

+  f^cosp  +  fufw)ds  +  puqtr  f*TA1L  R(R  sinp-x  cos  P)(-R  sin2  p 

+  x  sinp  cos  p  +  Rfu  cos  p  +  xfu  sin  p  +  fq  cos  p  +  fufq)ds 

+  pvviT  f^rAlL  R(R  sin  p  -  x  cos  P)h  wds 


+  pq*  /0TAIL  R(R  sinp-x  cosp)hqds  (26) 


\ 
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Since  the  only  potential  axial  force  of  practical  significance  is  the  axial  apparent 
mass,  ail  the  terms  in  p  Q,  except  the  term  containing  u,  were  dropped  when  de¬ 
ter  ming  Fa  . 


HYDRODYNAMIC  COEFFICIENTS 


The  expressions  for  FT,  FA,  and  M  in  terms  of  the  SNAME  (see  Ref.  7) 
hydrodynamic  nomenclature  are 

Ft  =  Z^-i-pi’w+Z^ip^q+Z^lp^uw+Z^ip^uq 

Fa  (apparent  mass)  =  X^pi3u 

M  =  ipf4w  +  M^ipf5q  +  M^ipf3uw  +  M^pf4uq 


(27) 


Equating  the  expressions  for  FT,  FA,  and  M  given  by  Eq.  26  and  27,  and  solving 
for  the  hydrodynamic  coefficients,  the  following  is  obtained: 


Zw  = 


z«  - 


it  J^TAIL  R  cos  phw  ds 
1/3 

TT  fJ™L  R  cos  Ph  ds 


1/4 

2* 


tt  J^TA!L  R  cos  (3(-  sin  p  cos  P  -  fu  sin  p  +  f  w  cos  P  +  fufw)da 


Z  = 

w 


it2 


K  - 


tt  j^AIL  R  cos  (3(-R  sin2  (3  +  x  sin  p  cos  P  +  Rf  u  cos  (3  +  xf  u  sin  p  +  f q  cos  (3  +  fufq)ds 


I«3 


(  tt  R(R  sinp  -  x  cos(3)hwds 


M  -  = 

W 


m'.  = 

4 


if4 


ttJ^tailR(R  sin  p  -  x  cos  P)hqds 


1*5 


M'  = 


tt  /qSTAIL  R(R  sin  p  -  x  cos  p)(-  sin  p  cos  p  -  fu  sin  p  +  fw  cos  p  +  fufw  )ds 


if3 


12 


NAVWEPS  REPORT  8  799 


it  fnTAIL  R(R  sin(3-x  cos  (3 )(- R  sin2  (3+x  sin  (3  cos  {3+Rfu  cos  (3+xfu  sin  (3+fa  cos  (3+f  f  )ds 

Ml  * 0  M  u  q 


1#4 


x'.  = 

u 


2tt  /JTAIL  R  sin  (3hu  ds 


i*3 


(28) 


EVALUATION  OF  INTEGRALS 


The  integrals  in  Eq.  28  that  express  the  hydrodynamic  coefficients  are  functions 
of  the  body  geometry  and  the  functions  fu(s),  fw(s),  fq(s),  hu(s),  hw(s),  and  hq(s). 
From  Eq.  8,  9*  and  12  it  is  seen  that 


(v  )  =  uf  (s) 

us  n  -O  u 

<Vws  >nO  =  wfw<S>  COS  0 
( Vqs  U=  COS0 

(*  )  =  u(F  )  ^  =  uh  (s) 

'  u  'n=Q  '  u  n^O  u 

(4>w  )n  0  =  w(Fw>n=0  COS  0  =  whw<S)  COS  9 


<<t>q)n^)=  q(F  )BSOCO8  0  =  qhq(s)  cose 

Hence  the  functions  f  u ,  fw,  and  fq  describe  the  variation  of  the  fluid- velocity 
component,  in  the  s  direction  at  the  body  surface  and  in  the  meridian  plane  0=0, 
for  unit  values  of  u,  w,  and  q.  Further,  the  functions  hu,  hw,  and  hq  describe 
the  variation  of  the  potential  at  the  body  surface  in  the  meridian  plane  0=0  for 
unit  values  of  u,  w,  and  q.  As  discussed  earlier,  these  functions  are  given  by  the 
extended  Douglas  program.  They  are  obtained  from  the  present3  version  of  this 
program  in  the  following  manner.  To  obtain  fu  and  hu,  the  following  items  are 
checked  on  the  control  cards:  (1)  surface  of  revolution,  (2)  perturbations  only, 
and  (3)  potential  computed.  In  the  outputs,  T  =  fu  and  PHI  =  hu .  To  obtain  f*  and 
hw,  the  items  checked  on  the  control  cards  are  (1)  crossflow,  (2)  perturbations 
only,  and  (3)  potential  computed.  In  the  outputs,  T2  =  -fw  and  PHI  =  -hw.  To  ob¬ 
tain  fq  and  hq,  the  items  checked  on  the  control  cards  are  (1)  nonuniform  onset 
flow  (computed),  (2)  perturbations  only,  and  (3)  potential  computed.  In  the  out¬ 
puts,  T2  =  f  and  PHI  =  hq. 


3A  new  program  for  calculating  the  potential  flow  nbout  axisymmetric  bodies  is  currently  being  developed  at  the 
Douglas  Aircraft  Co.  This  program  is  written  in  FORTRAN  IV,  in  contrast  to  the  present  machine -language  version 
that  operates  in  the  SAMSON  system. 
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A  computer  program  was  written  at  this  Station  to  perform  the  in¬ 
tegrations  indicated  in  Eq.  28.  Since  the  values  of  fu,  fw  ,  f q  and  hu, 
hw  ,  hq  obtained  from  the  Douglas  program  will  not  normally  occur  at 
equally  spaced  values  of  s,  values  of  the  integrands  are  calculated  at 
the  points  where  fu,  f w  ,  fq  and  hu,  hw,  hq  are  known  from  the  Douglas 
program  solution.  Values  of  the  integrands  at  equally  spaced  values  of 
s  are  then  determined  by  three-point  Lagrange  interpolation,  and  the 
integrals  art  evaluated  using  Simpson's  rule.  An  increasingly  larger 
number  of  points  is  used  in  the  Simpson- rule  integration,  until  the 
values  of  the  hydrodynamic  coefficients  reach  an  essentially  constant 
value. 


A  COMPARISON  WITH  EXPERIMENT 

Since  the  analysis  thus  far  has  been  limited  to  a  strictly  potential 
flow  solution,  a  comparison  with  experimental  data  from  a  blunt-based 
body  is  the  only  valid  one  that  can  be  made,  because  the  boundary 
layers  on  such  bodies  are  usually  so  thin  that  they  may  be  neglected. 
(For  a  discussion  of  this  area,  see  the  Introduction.  )  Unfortunately, 
the  only  blunt-based  bodies  for  which  complete  sets  of  measured 
hydrodynamic  coefficients  exist  are  the  Polaris  missiles.  The  A1 
configuration  was  chosen  for  the  comparison  be  :ause  it  was  hoped  that 
the  hydrodynamic  coefficients  for  this  early  version  of  the  missile 
could  be  declassified.  This,  however,  proved  impossible.  Therefore, 
the  comparison  is  made  on  the  basis  of  the  ratio  of  the  coefficient  de¬ 
termined  from  the  theoretical  calculations  to  the  measured  coefficient. 

In  calculating  the  hydrodynamic  coefficients,  it  was  necessary  to 
account  for  the  fact  that  on  a  blunt-  or  flat-based  body  the  flow  sepa¬ 
rates  from  the  body  at  the  base.  For  an  observer  traveling  with  an 
axisymmetric  blunt-based  body  having  (1)  a  large  forward  velocity,  (2) 
a  small  angle  of  attack,  and  (3)  small  pitch  rate,  the  flow  behind  the 
body  travels  rearward  almost  as  though  the  body  extended  back  for 
for  some  distance  behind  its  true  termination.  Therefore,  in  comput¬ 
ing  the  values  of  fu,  fw,  and  f q ,  a  simple  model  which  assumes  the 
body  to  extend  back  indefinitely  at  its  base  diameter  was  used.  Of 
course,  in  calculating  the  values  of  Z^»  Zq,  Mw  *  and  Mq  from  these 
values  of  fu,  fw,  and  fq,  the  integrations  were  performed  over  the  true 
surface  of  the  body  only.  In  considering  a  model  to  use  for  calculat¬ 
ing  the  values  of  hw  and  hq,  it  appeared  likely  that  the  pressure  field 
created  in  the  fluid  due  to  a  transverse  or  angular  acceleration  of  the 
body  would  fairly  strongly  sense  the  presence  of  the  abrupt  termination 
of  the  body  at  the  base.  Hence  the  values  of  hw  and  hq  used  in  calcu¬ 
lating  the  values  of  Z%-v>  Z<j  ,  M^»  Mq  were  obtained  with  a  model  which 
assumes  that  the  body  ends  at  its  true  base  termination,  but  does  not 
have  a  solid  surface  across  the  base.  Although  several  other  models 
were  tried  for  calculating  values  of  fu,  fvv,  and  fq  and  values  of  hw  and 
hq ,  the  two  models  discussed  above,  which  are  the  simplest  possible 
from  a  computational  standpoint,  were  found  to  give  the  best  results. 
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The  question  arises  whether  these  representations  of  the  base  flow, 
which  worked  well  for  the  A1  Polaris  missile  configuration,  will  work 
well  for  bodies  of  greatly  different  shape.  Owing  to  a  lack  of  experi¬ 
mental  data,  this  question  cannot  be  answered  at  this  time  and  will  have 
to  await  further  experimental  studies. 

A  comparison  of  the  theoretical  and  measured  hydrodynamic  coeffi¬ 
cients  of  the  A1  Polaris  missile  is  given  in  Table  1. 


TABLE  1.  Ratios  of  Theoretical  (T)  to 
Measured  (M)  Values  of  Hydrodynamic 
Coefficients,  A!  Polaris  Missile 


Coefficient 

Numerical 

value  of  ratio 

Error  in 

theoretical 

coefficient 

(K*h 

<Zw>M 

0.932 

-6.8% 

<z;>r 

(Zi>M 

0.927 

-7.3% 

Kvh- 

(<>m 

0.929 

-7.1% 

<Mq>T 

<Mq>M 

1.047 

+4.7% 

<ZJv>T 

<Zw)m 

1.023 

+2.3% 

<z4>r 

<z4>m 

0.935 

-6.5% 

(mW>T 

(Mw)m 

0.992 

-0.8% 

Kjh- 

(m4>m 

0.961 

-3.9% 
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Values  given  by  the  new  theoretical  approach  presented  in  this  report 
are  compared  with  the  measured4  values  given  in  Ref.  8.  As  mentioned 
above,  the  comparison  is  made  on  the  basis  of  the  ratio  of  theoretical 
to  measured  values. 

The  largest  error  is  7.3%.  Considering  the  probable  experimental 
error  in  the  measured  hydrodynamic  coefficients,  this  new  theoretical 
approach  appears  to  give  values  for  blunt-based  bodies  that  are  ade¬ 
quate  for  most  engineering  work.  Naturally,  this  supposition  is  tenta¬ 
tive,  since  the  comparison  is  limited  to  only  one  body  shape. 

A  comparison  similar  to  that  shown  in  Table  1  w  s  made  using 
slender  body  theory.  In  Table  1,  the  mean  of  the  absolute  values  of  the 
percent  error  is  4.9;  for  slender  body  theory,  it  is  11.7.  In  Table  1, 
the  maximum  percent  error  is  7.3;  for  slender  body  theory,  it  is  23.5. 

It  thus  appears  that  the  new  theoretical  approach  yields  a  substantial 
improvement  over  slender  body  theory. 


DISCUSSION  AND  CONCLUSIONS 

V  new  method  has  been  presented  for  calculating  the  hydrodynamic, 
coefficients  of  axisymmetric  bodies.  This  method  utilizes  an  extended 
version  of  the  Douglas  potential  flow  computer  program  for  axisymmetric 
bodies  to  determine  those  velocities  and  potentials  at  the  body  surface 
which  are  required  to  solve  for  the  forces  and  moment  acting  on  the 
body.  The  development  here  has  been  restricted  to  blunt-  or  flat-based 
bodies  where  the  boundary  layer  on  the  body  may  be  neglected.  It  is 
P1  anned  to  extend  this  work  to  bodies  having  streamlined  aftersection 
shapes  where  the  thick  boundary  layer  must  be  accounted  for. 

Based  on  a  comparison  with  experiments  limited  to  one  body  shape, 
it  may  tentatively  be  concluded  that  the  new  method  yields  hydrodynamic 
coefficient  values  for  blunt-based  bodies  that  are  adequate  for  most  en¬ 
gineering  purposes.  The  successful  application  to  a  blunt-based  body 
indicates  that  the  method  might  eventually  be  extended  to  streamlined 
bodies. 


^The  Zyy  and  Mw  coefficients  were  measured  in  the  water  tunnel  of  the  Ordnance  Research  Labora¬ 
tory,  Pennsylvania  State  University.  The  remaining  coefficients  were  measured  on  the  pitch-and-heave 
oscillator  at  the  David  Taylor  Model  Basin,  Washington,  D.C. 
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